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Controlling the spin orientation of photoexcited electrons by symmetry breaking 

Lan Qingjj Yang Song, and Hanan Derjjj 
Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627 

We study reflection of optically spin-oriented hot electrons as a means to probe the semiconductor 
crystal symmetry and its intimate relation with the spin-orbit coupling. The symmetry breaking by 
reflection manifests itself by tipping the net-spin vector of the photoexcited electrons out of the light 
propagation direction. The tipping angle and the pointing direction of the net-spin vector are set by 
the crystal-induced spin precession, momentum alignment and spin-momentum correlation of the 
initial photoexcited electron population. We examine non-magnetic semiconductor heterostructures 
and semiconductor/ferromagnet systems and show the unique signatures of these effects. 

PACS numbers: 72.25.Fe, 72.25. Mk, 72.25.Rb, 71.70.Ej 



The intimate relation between the spin-orbit coupling 
of crystals and their symmetry is a theme of on-going 
research for more than a half-century [THHj. The rapid 
progress of spintronics research has provided new pow- 
erful techniques to study this relation [TS]. In semicon- 
ductors it is readily seen in the valence band energy dis- 
persion [11 12], in spin relaxation of electrons [3], or in 
optical selection rules |4] [5] . In magnetic materials this 
relation sets the magnetocrystalline anisotropy and mag- 
netostriction constants [TSIIIH. In this letter, we study 
intriguing results of this relation which depend on the 
spin-momentum correlation and coherent spin precession 
of photoexcited electrons in bulk semiconductors. 

Using symmetry breaking by reflection we show that 
the net-spin vector of photoexcited electrons is tipped 
away from the light propagation axis. We explain the 
dependence of the tipping on spin-momentum correla- 
tion and coherent spin precession, and then calculate the 
tipping angles with partial, complete, nonmagnetic and 
magnetic reflections. In prospective experiments, one can 
measure signatures of the correlation and coherent pre- 
cession by probing the tipped net-spin vector rather than 
the coherent distribution of photoexcited electrons. One 
clear advantage is that the net-spin vector decays orders 
of magnitude slower than the coherence time. 

We first explain the correlation tipping phenomenon 
in a photoexcited direct gap semiconductor. Figure ma) 
shows the momentum alignment and spin-momentum 
correlation of electrons immediately after photoexcita- 
tion. The shown distributions are compactly derived by 
using the spherical model at the top of the valence band 
[T] and s-type states in the conduction band. The density 
matrix of photoexcited electrons is then [U [S] , 

V{t, k) = |l + ia(t) [3|e • k|^- 1]! + 2S{t){& ■ p) 

+ i/3(i)[3(6-.k)(p.k)-(<^.p)]}F(i,fc). (1) 

k (e) is the unit vector in the direction of electron mo- 
mentum (light polarization). & (Z) denotes the Pauli 
matrix vector (the 2x2 unit matrix). The photon an- 
gular momentum unit vector is defined by p = i e x e* . 




FIG. 1: (a) Momentum and spin distributions of photoex- 
cited electrons following transitions with (i) heavy holes and 
(ii) light holes (excitation with heavy holes is stronger due 
to their larger density of states). The length of a line from 
the center represents the relative population of excited elec- 
trons with momentum along the arrow's direction. The thick 
arrows on the edge line represent the correlated spins, (b) 
Correlation induced tipping by partial reflection. The net- 
spin vector is tipped away form p due to missing transmitted 
electrons with spins mostly along the interface normal n. (c) 
Spin precession of electrons immediately after photoexcita- 
tion. Electrons moving in opposite directions have the same 
initial spin direction but precess at opposite angles, (d) Pre- 
cession induced tipping by complete reflection. The net-spin 
vector is tipped away from p due to rephasing (see text). 



For linearly polarized light p = and for circularly po- 
larized light, |e • kl = [1 — (p • k)^]/2. The parame- 
ters a, S and /3 are, respectively, measures of momentum 
alignment, of average spin and of spin- momentum corre- 
lation. S(t) decays exponentially with the spin relaxation 
time from initial value Sq ~ —1/4. a(t) and /3(t) decay 



exponentially with the much shorter momentum relax- 
ation time from their initial values ao ~ /3o ~ il |16j . 
The lower/upper sign is for transitions with heavy/light 
holes. Finally, F{t, k) relates to the density of excited 
electrons. Figure [l|b) elucidates the correlation induced 
tipping effect for electrons that are generated with heavy 
holes. This effect is governed by transmission of electrons 
immediately after excitation (prior to the decay of /3{t) 
in Eq. (fl])). For example, consider electrons that prop- 
agate parallel and perpendicular to the reflection plane 
(shown in Fig. W[h) by the two spins on the edges of the 
dash line and the two spins along ±n). Their net spin is 
not coUinear with p due to missing electrons along —n. 
Similar explanation holds for other spin and momentum 
directions. The overall effect is that the net-spin vector 
changes from 5oP to Sq [(1 — Sa)p — Sh{p ■ n)n] where Sa 
and Sb are measures of the transmission amplitude. 

The second tipping phenomenon results from intrinsic 
spin precession of electrons in semiconductors without in- 
version symmetry [3]. The precession is due to a torque 
exerted by the electron's effective magnetic field whose 
components are Bi oc ki(k'^ — fc^) where {ki, kj, km} are 
the electron's wavevector components along the crystal- 
lographic axes. Figure llTc) shows the intrinsic spin pre- 
cession of photoexcited electrons immediately after gen- 
eration. Electrons that move in opposite directions have 
similar initial spin direction but precess at opposite an- 
gles. The net angle precession of the pair is zero on av- 
erage. This picture is changed by specular reflection of 
one of the electrons off an interface as shown in Fig. flld) . 
Here, kz and Bz of an electron flip signs after specular re- 
flection whereas its spin components {S^, Sy, Sz} remain 
unchanged. This causes Sx and Sy to 're-phase' to some 
degree due to the flip of Bz whereas Sz keeps its pre- 
cession unperturbedly. The net effect is that transverse 
spin components with respect to the normal of the reflec- 
tion plane decay slower than the longitudinal component. 
This effect is robust in both complete or partial reflec- 
tion and it can be further amplified by multiple interface 
reflections. Scattering events reduce the magnitude of 
the net-spin vector but they do not change its direction. 
We will show that in spite of a notable Dyakonov-Perel 
spin relaxation of hot electrons the net-spin vector after 
energy relaxation remains measurable. 

The time evolution of photoexcited electrons is stud- 
ied by extensive Monte Carlo simulations [T7] . The initial 
wavevector and spin directions are randomized according 
to the distributions in Eq. (IT]). The initial electron den- 
sity follows the light attenuation profile. We use the ef- 
fective mass approximation in calculating electron trans- 
port and quantum mechanical transmission probabilities 
across interfaces. Effective masses, band-gaps and band 
offsets are taken from Ref. [18]. Momentum and energy 
relaxations of hot electrons are simulated by the Frohlich 
interaction [TS] . Between scattering and reflection events 
the spin precesses about its intrinsic field. The simulation 
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FIG. 2: (a) Simulated excitation geometry and heterostruc- 
ture setups (i) and (ii). n is along the z crystal axis, p is in 
the x-z crystal plane (cosS = p • n). Sr is the remained net- 
spin vector in the GaAs region after energy relaxation. A^ is 
the tipping angle of S^ away from the optically injected spin 
direction, (b) Ar vs Q for various setups and photon energies, 
(c) Ar (left) and Sr (right) vs photon energy in setup (i). The 
decay of Sr is due to hot-electron spin relaxation, (d) Spin 
evolution along crystal axes in the GaAs region in setup (i). 
(e) Ar vs photon energy in setup (ii) with/without alignment 
and correlation, (c)-(e) correspond to 6 = 7r/4. 



ends when electrons reach the bottom of the conduction 
band (typically after f ps) . Spin relaxation at later times 
occurs on >1 ns time scales [12 . 

We first study the tipping effects in non-magnetic het- 
erostructures. To distinguish between correlation and 
precession induced tipping we employ two setups of the 
heterostructure AlAs/GaAs(50 nm)/AlxGai_xAs(50 nm) 
/Alo.osGao.gsAs [2^. The structure, shown in Fig. [2ja), 
ensures that confinement effects are negligible and all re- 
gions are bulk in nature. Setup (i) includes a high in- 
ner barrier (x=0.25) for which the potential steps are 
Ex = 230 meV h E2 = 185 meV (see Fig. ^a)). The 
tipping angle is governed by spin precession of hot elec- 
trons in the GaAs region prior to thermalization while 
bouncing back and forth from the potential walls. Setup 
(ii) includes a shallower inner barrier (x=O.I) for which 
the potential steps are E\ = 90 meV Sz E2 = 4:5 meV. For 
certain photon energies there is a favored net transmis- 
sion from the GaAs region to the Alo.05Gao.95As region. 
The spins of transmitted electrons are mostly aligned 
with the interface normal rather than the optically in- 
jected spin direction. In both setups, quantum tunneling 
across the 50 nm inner barrier is negligible. 

Figure [2]jb) shows the dependence of the tipping angle 
on the light propagation direction for various photon en- 
ergies. Ar is the tipping angle of the net-spin vector in 
the GaAs region after energy relaxation. 9 is the angle 
between the photon angular momentum (p) and the in- 
terface normal (n) as shown in Fig. [2ja) . We first focus 



on the results of setup (i). The tipping angle is max- 
imized when light propagates along the (101) & (Oil) 
crystallographic axes {9 = 7r/4) due to the faster spin 
precession along these directions. Similarly, the tipping 
angles increase with photon energies due to the enhanced 
precession of electrons. Figure [2]^c) shows the tipping an- 
gle as a function of photon energy when 9 = 7r/4. The 
net-spin vector after energy relaxation remains measur- 
able in all of the studied excitation range. Figure W[d) 
shows the time evolution of the net-spin along the crystal- 
lographic axes for 9 = tt/4 and photon energy of 1.75 eV. 
The aforementioned rephasing effect is seen in the slower 
decay of Sx during the energy relaxation in the first 0.5 ps 
(the z component of the intrinsic field changes direction 
with each reflection). At later times (>1 ps) the preces- 
sion frequency and spin relaxation are much slower. 

Figure pie) shows correlation induced tipping angles 
as a function of photon energy when 9 = 7r/4. As 
a reference, we also simulate a case where momentum 
alignment and spin-momentum correlation are neglected 
(dash line). In the latter case, only spin precession dur- 
ing the initial coherent phase induces the tipping angle. 
The first sharp increase is reached when electrons have 
enough energy to cross to the barrier region. The miss- 
ing transmitted electrons have spin components mostly 
along n (generated with heavy holes). This increase is 
then suppressed (plateau region) by the backward trans- 
mission from the Alo.05Gao.95As [21]. The following drop 
is due to the transmitted hot electrons generated with 
light holes from both sides of the barrier (see Fig. W[a) 
for the differences in spin-momentum correlation of the 
heavy and light holes cases). Further increase of photon 
energy suppresses the correlation induced tipping mech- 
anism due to the increased precession rates. 

Partial reflection off a ferromagnet provides unique sig- 
natures of momentum alignment and spin-momentum 
correlation that are interwoven with the magnetization 
direction of the ferromagnet. We perform detailed Monte 
Carlo simulations for the heterostructure in Fig. IsFa). 
It consists of Fe/GaAs(150nm)/Alo.3Gao.7As(50nm) 
/Alo.25Gao.75As(50nm)/AlAs. For certain photon ener- 
gies reflected hot electrons from the Fe/GaAs interface 
reach region I (Alo.25Gao.75As) only if they do not expe- 
rience energy relaxation in the GaAs region. In addition, 
the spins of these electrons do not precess in their short 
passage since their motion is along a crystallographic axis 
(n = z). Partial and spin selective reflection across the 
Fe/GaAs interface is modeled by a 0.5 eV high and 6 nm 
wide parabolic Schottky barrier [17] . 

Figure p[b) shows the tipping angle after energy re- 
laxation in region I as a function of the photon energy. 
At photon energies just below the band-gap of region I, 
transmission from the GaAs into region I across the inner 
barrier is possible only for electrons that are generated 
with heavy-holes and that are directed along n. This is 
shown in the low energy side of Fig. p[b). Since the spin 
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FIG. 3: (a) Simulated excitation geometry and heterostruc- 
ture. rh, n and mx n are along crystallographic axes (in-plane 
magnetization), (b) Tipping angle of the relaxed net-spin vec- 
tor in region I vs photon energy when the light propagation 
axis is 45° from fi. The contribution of alignment and cor- 
relation are clearly seen, (c) Relaxed net-spin along m x n 
in region I versus light propagation direction (the photon en- 
ergy is 1.85 eV). The polar (6) and azimuthal {(jj) angles are 
measured from n and ih x n, respectively. 



direction of these electrons is parallel to their wavevec- 
tor, the tipping angle is simply the angle between the 
optically injected direction and n. The step occurs when 
electrons are photoexcited in region I and the net-spin 
becomes aligned with the optically injected spin direc- 
tion (the transmitted electrons from the GaAs region are 
shadowed by the density of photoexcited electrons in re- 
gion I) . Figure pTc) shows the relaxed spin component 
along rh X fi in region I as a function of the light prop- 
agation direction for photon energy of 1.85 eV. This rel- 
atively small component is an entwined signature of the 
ferromagnet and the spin-momentum correlation. 

The simulated results in Figs. [2] & [3] are consistent 
with a simple quantitative analysis in which the pre- 
cession and correlation induced tipping are treated sep- 
arately. For both mechanisms, the average spin af- 
ter relaxation in the non-magnetic structure reads S = 
So [(1 — Sa)p — 6b{p ■ n)n]. In the optimal case (p • fi = 
l/\/2) the tipping angle is A « 5i,/2 where for the pre- 
cession induced tipping (setup (i) in Fig. [2]ja)) [T7] . 
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2h?E„ ^ 

3 n=0 



7p,„(/ia; ~ Eg - nELof 



(2) 



7p,ra=o ~ 0.06 and 7p,n>o ~ 0.08 are integration param- 
eters. Eg is the energy gap, and Uc is a measure of the 
spin-orbit coupling strength in the conduction band [3]. 
Elo is the longitudinal-optical phonon energy and t^q 
is the associated scattering time. N = [(fiw — Eg)/ELo\ 
denotes the number of phonon emissions prior to relax- 
ation to the bottom of the band. For the correlation 
induced tipping (setup (ii)) we get. 
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in the low photon energy regime (before the plateau in 
Fig. [2Fe)). 7c « 0.4 is an integration parameter and 



/3o corresponds to the correlation parameter of electrons 
generated from heavy holes, i? = [mhh{f>^ — Eg)]/[{me + 
mhh)Ei] — 1 where ruhh ("^e) is the heavy-hole (electron) 
effective mass and Ei is the barrier height (see Fig.l2ja)). 
Using Eqs. ([2|-([3]), the spin orbit coupling parameters 
ttc and /3o can be extracted from experiments |17j . In 
the magnetic case of Fig. Wla) we can also work out the 
tipping angles since precession effects are small (electrons 
reach the interface prior to momentum scattering). Here, 
the integrated surface density of the net-spin vector after 
magnetic reflection reads, 



Sr = no^S'op-no{(4AoS'o 



S^)P ~ 3ip ■ n)6^ n 



+ [Ai - A2 + i(l - 3|e • h\){S'l - 5^)]m 
-(4A3S'o + (5f)mxp + 3(p-n)(5f riix n}. (4) 

£~^ is the light absorption coefficient and no is the pho- 
toexcited electron density. Xi (Si) are isotropic and un- 
correlated (momentum aligned and spin correlated) in- 
tegration parameters that depend on the barrier trans- 
mission. The total reflection is spin independent and 
governed by Aq & Sq terms, the spin selective reflection 
by (Ai — A2) & {S1—S2) terms, and the magnetization in- 
duced torque by A3 & 6^ terms. The a//3 superscript in- 
dicates signatures of alignment/correlation. For detailed 
expressions of all terms see [T7]. Previous experimen- 
tal investigations of the ferromagnetic proximity effect 
[22t [23] and their ensuing theories [24H26] were focused 
on ferromagnetic signatures while ignoring the alignment 
and correlation of photoexcited electrons. Ciuti et al. 
have derived a reduced form of Eq. m in which 5f— (5,f = 
[21] ■ Our analysis shows that in properly designed struc- 
tures, correlation induced signatures are experimentally 
resolvable by the p • n amplitude dependence of the non- 
magnetic (n) and magnetic (m x n) spin components. 

In conclusion, reflections off non-magnetic semiconduc- 
tor heteroj unctions and semiconductor/ferromagnet in- 
terfaces have been shown to be a powerful tool to study 
coherent effects of the crystal symmetry and spin-orbit 
coupling. Tipping the net-spin vector out of the optically 
injected direction is a measure of these effects. The pre- 
dicted tipping angles are noticeable and can be probed, 
for example, via the photoluminescence of energy relaxed 
free excitons. The tipping angle corresponds to the an- 
gle at which the detected circular polarization is maxi- 
mal. Tunable parameters are the photon energy and light 
propagation direction. The tipped net-spin vector would 
ultimately evolve in hyperfine interaction and polarize 
the nuclear spin system. In this case, ultrafast decaying 
coherent effects that result from the crystal symmetry 
and spin-orbit coupling can be inferred by the 10^*^ slower 
dynamic polarization of the nuclear system. 

This work is supported by AFOSR and NSF (Contracts 
No. FA 9550-09-1-0493 and ECCS-0824075). We deeply 
thank Dr. Scott Crooker for helpful insights. 



* Electronic address: llan.qing@ rochester.edu| 
^ Also at Department of Electrical and Computer Engi- 
neering, University of Rochester, Rochester, NY 14627. 
[1] J. M. Luttinger and W. Kohn, Phys. Rev. 97, 869 (1955). 
[2] G. Dresselhaus, Phys. Rev. 100, 580 (1955). 
[3] M. I. Dyakonov, V. I. Perel, Sov. Phys. JETP 33, 1053 

(1971); Sov. Phys. Solid State 13, 3023 (1972). 
[4] V. D. Dymnikov, M. I. Dyakonov and N. I. Perel, Sov. 

Phys. JETP 44, 1252 (1976). 
[5] Optical Orientation, edited by F. Meier and B. P. Za- 

kharchenya, (North-Holland, New York, 1984). 
[6] P. Li and H. Dery, Phys. Rev. Lett. 105, 037204 (2010). 
[7] V. I. Zemskii, B. P. Zakharchenya, and D. N. Mirlin, 

JETP Lett. 24, 82 (1976). 
[8] V. L. Alperovich, V. L Behnicher, V. N. Novikov, and A. 

S. Terekhov, JETP Lett. 31, 546 (1980). 
[9] Y. A. Bychkov and E. I. Rashba, J. Phys. C 17, 6039 
(1984). 

[10] R. D. R. Bhat and J. E. Sipe, Phys. Rev. Lett. 85, 5432 
(2000). 

[11] J. Wunderlich, B. Kaestner, J. Sinova and T. Jungwirth, 
Phys. Rev. Lett. 94, 047204 (2005). 

[12] S. A. Crooker and D. L. Smith, Phys. Rev. Lett. 94, 
236601 (2005). 

[13] J. Kanamori, in Magnetism, edited by C. T. Rado and H. 
Suhl, (Academic Press, New York, 1963), Vol. I, p. 127. 

[14] T. Dietl, H. Ohno, and F. Matsukura, Phys. Rev. B 63, 
195205 (2001). 

[15] I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

[16] These values hold for photon energies below the energy 
gap between the conduction and split-off valence bands. 

[17] See supplementary material for simulation procedures, 
derivations of Eqs. Q-Q, and possible experiments. It 
also provides details that support our calculations. 

[18] Properties of Aluminium Gallium Arsenide, edited by S. 
Adachi, (INSPEC, London, 1993). 

[19] P. Y. Yu and M. Cardona, Fundamentals of Semiconduc- 
tors (Springer, Berlin, 2005), Ch. 5. 

[20] We choose Alo.05Gao.95As in the right region in order 
to enhance the net transmission out of the GaAs region. 
This choice also allows one to distinguish the free exciton 
luminescence of the GaAs region. 

[21] Application of an external electric field or the use of 
structures with built-in electric field (e.g., by continous 
variation of the Al content in the barrier region) would 
enhance the transmission from the GaAs region and it 
will also suppress the plateau appearance. 

[22] R. K. Kawakami Y. Kato, M. Hanson, I. Malajovich, J. 
M. Stephens, E. Johnston-Halperin, G. Salis, A. C. Cos- 
sard, D. D. Awschalom, Science 294, 131 (2001). 

[23] R. J. Epstein, I. Malajovich, R. K. Kawakami, Y. Chye, 
M. Hanson, P. M. Petroff, A. C. Gossard, and D. D. 
Awschalom, Phys. Rev. B 65, 121202 (2002). 

[24] C. Ciuti, J. P. McCuire, and L. J. Sham, Phys. Rev. Lett. 
89, 156601 (2002). 

[25] G. E. Bauer, A. Brataas, Y. Tserkovnyak, B. I. Halperin, 
M. Zwierzycki, and P. J. Kelly, Phys. Rev. Lett. 92, 
126601 (2004). 

[26] V. N. Cridnev, JETP Lett. 77, 187 (2003). 



